A complicated DufRng oscillator in the surface-electrode ion trap 
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The oscillation coupling and different nonlinear effects are observed in a single trapped ''"Ca"'" 
ion confined in our home-built surface-electrode trap (SET). The coupling and the nonlinearity 
are originated from the high-order multipole potentials due to different layouts and the fabrication 
asymmetry of the SET. We solve a complicated Dufhng equation with coupled oscillation terms 
by the multiple scale method, which fits the experimental values very well. Our investigation in 
the SET helps for exploring nonlinearity using currently available techniques and for suppressing 
instability of qubits in quantum information processing with trapped ions. 

PACS numbers: 05.45.-a, 37.10.Vz, 37.10.Ty 



I. INTRODUCTION 

The Duffing oscillator is generally used to describe non- 
linear dynamics in oscillating systems The corre- 
sponding Duffing equation models a damping and driven 
oscillator with more complicated behavior than simple 
harmonic motion, which can be used to exhibit chaos in 
dynamics and hysteresis in resonance 

On the other hand, the motion of trapped ions is highly 
controllable and can be employed to transfer quantum in- 
formation when cooled down to ground state [T]| . Since 
it is effectively approximated to be harmonic, the ion 
motion in a quadruple electromagnetic trap |12| can be 
regarded as a good mechanical oscillator, which may ex- 
hibit nonlinearity when driven to the nonlinear field. For 
example. Duffing nonlinear dynamics has been investi- 
gated in a single ion confined in the linear ion trap [l3| . 
The trap nonlinearity introduces instability in the motion 
of the ion, which should be avoided in most times, but 
can also be used in resonance rejection and parameter 
detection in mass spectrometry [ij, [l^ . Recent research 
also showed the feasibility of phonon lasers based on the 
nonlinearity of a single trapped ion under laser irradia- 
tion [III. 

We focus in this work on the nonlinearity in a home- 
built surface-electrode trap (SET). The SET, with capa- 
bility to localize and transport trapped ions in different 
potential wells, is a promisin g se tup for large-scale quan- 
tum information processing [17|. In comparison to con- 
ventional linear Paul traps, however, the reduced size and 
asymmetry in SET lead to stronger high-order multipole 
potentials (T8l - [20j . which affect the stability of the ion 
trapping. To solve the problem we have to understand 
the source and the strength of the nonlinearity. Due to 
complexity resulted from the high-order multipole poten- 
tials, the nonlinear effect in the SET cannot be simply 
described by the Duffing oscillator as for the linear trap. 



but an inhomogeneous-coupled Duffing oscillator involv- 
ing quadratic and cubic nonlinearities. We observed 
the nonlinearity experimentally in our SET, and by the 
method of multiple scales we derived an inhomogeneous- 
coupled Duffing oscillator to fit the experimental values, 
which shows that both the nonlinearity and axial-radial 
coupling exist in the case of the frequency resonance (i.e., 
around the regime of driving detuning being zero) . More- 
over, we show in the non-coupling case different nonlinear 
effects in different dimensions, which is due to different 
asymmetry in fabrication of the SET. 



II. 



EXPERIMENTAL SETUP AND IMAGES OF 
ION MOTION 
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FIG. 1: The layout of our home-buih SET with the top-right 
inset for the relevant energy levels of the trapped ''"Ca^ ion. 
The ''"Ca^ ion is cooled by the laser beam at 397 nm on the 
43^5*1/2 < — ^ transition, helped by another laser light 

at 866 nm as repumping. Here, the trapped ion is above the 
electrodes by 0.8 mm (See the bottom-right inset). 

Our home-built SET is a 500 //m-scaled planar trap 
with five electrodes for radial confinement, and fabricated 
by printed circuit board technology [21| . As shown in Fig. 
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1. the five electrodes eonsists of a eentral electrode, two 
radio-frequency (rf) electrodes and two outer segmented 
dc electrodes, where the rf electrodes, the central elec- 
trode and the gaps in between are of the same width of 
500 ^m. Each outer segmented electrode consists of five 
component electrodes, i.e., a middle electrode, two con- 
trol electrodes and two end electrodes. The widths of the 
control electrodes and end electrodes in the segment are 
1.5 mm and the middle electrode is 1 mm wide. The gap 
in the segmented electrode is of 500 /xm width. When 
the SET works, the trapped '*°Ca"'" ion stays above the 
electrodes by 0.8 mm, and the pseudopotential trapping 
depth is below 1 eV with rf amplitude J7(0-peak)= 400 
V and rf frequency Q,/2tt = 15 MHz. The voltage on 
the four end electrodes is Vec = 40 V but zero on other 
electrodes. 
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FIG. 2: Schematic of our experimental setup. The main com- 
ponents are explained in the text. 



The experimental setup is plotted in Fig. 2, where 
the ultraviolet radiation at 397 nm excites 45^5*1/2 < — ^ 
4p^Pi/2 transition by a grating stabilized laser diode with 
power up to 30 mW and linewidth less than 2 MHz. An- 
other grating stabilized laser diode at 866 nm with power 
up to 100 mW and linewidth less than 5 MHz excites 
the ?>(PD^/2 < — > 4p^Pi/2 transition. The frequencies of 
the both laser diodes have been calibrated to the wave- 
length meter (HighFiness, WS-7). Typical laser powers 
at the trap center are 50 /iW for 397 nm in red detun- 
ing (-80 MHz) and 500 for 866 nm in carrier tran- 
sition. In our SET, the single '*°Ca+ ion is laser cooled 
and stably confined, which is monitored by photon scat- 
tering collected by both an electron-multiplying charge- 
coupled-device (EMCCD) camera and a photomultiplier 
tube (PMT). Outside the vacuum chamber, the electrical 
connections immediately encounter a "filter box", which 
provides low-pass filtering of the voltages applied to the 
electrode. An additional drive force is electrically con- 
nected to one of the middle electrodes behind the filter 
box, which provides an excitation to drive the ion away 
from equilibrium. Due to the design of our SET system. 



the motion of the ion is detected only in the xz plane 
by the EMCCD. As a result, what we study throughout 
the work is the oscillation along the axial direction {z- 
axis) and the radial direction (x-axis), whose harmonic 
frequencies are, respectively, a;oz/27r = 191.7 kHz and 
a;o2;/27r ~ 425 kHz. Moreover, since the harmonic fre- 
quency in j/-axis is LOoy/^n = 925 kHz, much bigger than 
in other axes, the ion can be regarded as a very tight con- 
finement in y direction. We have suppressed the micro- 
motion by the rf-photon cross correlation compensation 
j2^ , which yield cooling of the ion down to the temper- 
ature below 10 mK. 

To study the nonlinear mechanical response, we drive 
the ion to the nonlinear regime by a small oscillating volt- 
age, i.e., V = 7Y, applied to one of the middle electrodes. 
We slowly increase the driving frequency with the scan 
step 0.1 kHz, from 189.0 kHz to 433.0 kHz (the positive 
scan), the ion oscillates first along the z-axis and then 
turns to the x-axis for oscillation. The particularly in- 
teresting observation is the simultaneous responses, i.e., 
a rectangle trajectory, in both x- and z- axes when the 
sweep is close to the harmonic resonator frequency in ei- 
ther of the axes. Similar behavior is also found in the 
negative scan. 
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FIG. 3: Time-averaged images of a single trapped ion taken 
at different drive frequencies, with (a) the positive scan and 
(b) the negative scan. The prime oscillation is along the axial 
direction and the scan step is 0.1 kHz. The rectangles ap- 
pear in the images at resonance frequency due to axial-radial 
coupling. 

We measure the ion oscillation by taking time-averaged 
images from the EMCCD. In Fig. 3, seventeen such 
images for different drive frequencies are presented for 
positive and negative scans, respectively. For the ion 
originally oscillating in the z-axis, we slowly scan the 
drive frequency Wz across the harmonic resonance at tdoz ■ 
When the detuning a = — ujqz approaches zero, the 
rectangle trajectory appears, implying a coupled motion 
between x- and z-axes due to axial-radial coupling (ex- 
plained later). For a more clarified observation, we scan 
with smaller steps around the regime cr = 0, as shown in 
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FIG. 4: The coupled oscillation in the time-average images 
for the positive scan. Comparing to Fig. 3, the scan with 
smaller step (0.01 kHz) presents clearer pictures for the rect- 
angle trajectories from appearance to disappearance, where 
sixteen images are selected in series with the frequency step 
0.05 kHz. 



Fig. 4 which gives us an accurate range from the appear- 
ance of the rectangle to the disappearance. 



III. THEORETICAL MODEL 

To understand the observation above, we have to con- 
sider the multipole potential in the SET, which is given 
by HI 



9ijMjYj {x,y,z) , 



(1) 



where the subscript i labels different electrodes and the 
subscript j is for the spherical harmonics Yj(x,y,z). 
Both Yj (XjV, z) and the related parameters Mj are de- 
fined in [23]. gij is the weight factor for different elec- 
trodes. In this treatment, the initial equilibrium posi- 
tion of the single trapped ion is defined as the origin of 
the coordinates. The five-wire SET generally consists of 
quadrupole and hcxapolc potentials [20] . Considering the 
defect in our SET, we also involve octopole potential in 
our treatment. Following the definition in [23|, we have 
the subscripts j = 5, • • • ,25 where j from 5 to 9, from 
10 to 16 and from 17 to 25 correspond, respectively, to 
the quadrupole, hexapole and octopole potentials. For 
different potentials =Vi + Ui cos(rit) applied, respec- 
tively, to N electrodes, where Vi is the dc voltage on the 
electrode i and Ui cos{Qt) represents the rf voltage Ui on 
the electrode i driven at frequency 17, we rewrite Eq. (1) 
for the dc potential ^dc and the rf potential <i>r/ as 

<^ac {x, y, z) = E -is yj*M,Y, {x, y, z) , 
<^rf (x, y, z) = E -Is U* cosmM.Y, {x, y,z), ^ > 

where we have used U* = X^i^i Uigij and V* = 

Moreover, we have the ID motional equation for the 
trapped ion [23 ]. 



with ^ = x,y and z. Combining Eq. (2) with Eq. (3), 
we obtain the equation of motion in the z direction. 




+aizy + a^zx 



-I- a^zxy -) 



+ 032: 
ajzy - 



-I- q;2iz 2/ + a22Z X 
aszx = kz cosiwzt), 
(4) 

where x, y and z represent, respectively, the displacement 
of the ion from the equilibrium position in the three di- 
mensions, ^ is the linear damping parameter originated 
from the recoil due to photon absorption, kz is driving 
amplitude. The detailed expressions of the nonlinear co- 
efficients ai [i = 2,21,22,3,4,5,6,7,8) can be found in 
Appendix I. Compared to the Duffing oscillator in [isj . 
Eq. (4) is a complicated Duffing oscillator, containing 
additional coupled-motion terms. 

Using the method of multiple scales [2F 
steady-state solution to Eq. 



we obtain the 



3a^ 



("3 



(4) as 
-2asa2xc'^ 



3w2^c 



3a2 



2q!5c2 
3a2 



)± 



(5) 



where the nonlinear coefficients a2x and a2y are relevant 
to the coupled motion along x- and y-axes. c, b and a 
are the response amplitudes, respectively, in x, y and z 
directions, ujox/^tt and a;oy/27r represent the harmonic 
frequencies in x- and y-axes. For more clarification, we 
define a parameter atotai as 



atotai = "3 + Aa2 + Aa, 



(6) 



where originates from the cubic nonlinearity, Aa2 — 
^-^T^ represents the nonlinear coefficient that comes 



9"n 



-2Q8Ct2xC" 



from quadric nonlinearity, and Aa - 

-L + corresponds to the nonhnear 



dispersion relevant to the coupled motion. Substituting 
cttotai into Eq. (5), we obtain 



3a 



total 2 I 

a ± 



8uJoz 



(7) 



IV. DISCUSSION ABOUT THE 
NONLINEARITY AND COUPLING 

In our home-built SET, since the harmonic frequency 
in y-axis is much bigger than in other axes, the ion is con- 
fined very tightly in y direction, which leads to a reason- 
able assumption h/a <^1. As a result, the coupled term 
Aa is reduced to Aa = xcVa^ with x = ^^f^ + 
Moreover, as and Aa2 in Eq. (6) are nothing to do with 
the coupled motion and their sum + Aa2 can be mea- 
sured experimentally by 



8^02 (Tri 



3(a3 -I- Aa2) ' 



(8) 
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with the maximal amphtude Um, and the maximal detun- 
ing am in the non-coupling case. As a result, Eq. (7) 
is reduced to a steady-state solution to the amplitude of 
the response a with respect to the drive detuning a for 
the known driving force amplitude kz, 



3a^ 



■(as -I- Aa2 




4^0. a' 



(9) 
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FIG. 5: The measured and calculated amplitudes with re- 
spect to the driving detuning, where the black dashed lines 
indicate the detunings for maximum amplitudes of oscilla- 
tions, (a) The measured values correspond to the situation in 
Fig. (3). Calculation (black solid curve) by Eq. (9) without 
the coupling term fits most of the measured values with pos- 
itive scan (red stars) and negative scan (green crosses). The 
scan step is 0.1 kHz. (b) The measured values correspond to 
the situation in Fig. (4). Calculation (light blue solid line) 
by Eq. (9) with the coupling term fits the measured values 
with positive slow scan (red stars) around zero detuning. The 
scan step is 0.01 kHz, and the blue stars are for the measured 
values of the vibrational amplitude c in x-axis. As a compari- 
son, the calculation without the coupling term is plotted (the 
black solid curve). 

In our experiment, is measured via ion response in 
the linear regime, — 0.075 x 10^ Hz^m is obtained by 
observing the ion displacement versus the middle elec- 
trode vohage HI, -t- Aaz = 0.1959 x Hz 
is measured using the observed dependence of am on the 
maximal detuning tT,„. We evaluate \i = 177.1 Hz using 
the relation a„i = kz/ {2imjJoz). The comparison in the 
non-coupling case between the measured and calculated 
values of a and a is made in Fig. 5(a), where Eq. (9) 



without the coupling term Aa (the black solid curve) 
can fit most experimental values for both the positive 
and negative scans (red stars and green crosses, respec- 
tively). In this situation, the vibrational amplitude c in 
X-axis is negligible. Some experimental values around 
(7 = 0, which are not fitted well by the solid curve, are 
actually relevant to the case of coupled motion. To be 
more clarified, we scan the region around a ~ with 
smaller step than in Fig. 5(a). The fitting by consid- 
ering the coupling term Aa in our calculation can fully 
cover the measured data, as shown in Fig. 5(b). In such 
a case, we find that the vibrational amplitude c in a;-axis 
is visible, which is excited by the energy transfer from 
2-axis due to motional coupling. This energy transferred 
from z-axis to x-axis is nearly constant in the adiabatic 
operation so that we obtain ~ 4.5 x 10^® Hz^/m^. 
Fig. 5(b) also shows that the motional coupling stops 
when (t/2tt approaches 0.25 kHz. We see that a goes up 
to a maximum with c dropping to zero, implying that 
the system returns to the non-coupling case. Therefore, 
the vibrational trajectories imaged in Figs. 3 and 4 can 
be fully understood by the complicated Duffing oscillator 
with and without the term for coupled motion. 
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FIG. 6: The measured and calculated amplitudes in a;-axis 
with respect to the driving detuning, where the scan step 
is 1 kHz and the measurement is made in the no-coupling 
case. Calculation (black solid curve) fits most of the measured 
values with positive scan (red stars) and negative scan (green 
crosses) . 

Moreover, we also checked nonlinear effects in differ- 
ent directions in our SET by applying the drive on the 
X-axis and repeating the experimental steps as above for 
z-axis. It is physically evident that the behavior can be 
described by a slight modification of Eq. (4) by exchang- 
ing z and X, and replacing kz and uJz by k^ and uj^- As an 
example, we only present the non-coupling case in Fig. 
6, where the measured data is fitted well by the steady 
solution to the complicated Duffing oscillator with dif- 
ferent nonlinear coefficients and different damping rates 
from in Fig. 5(a). In comparison to the shape of the 
curve in Fig. 5(a), the oscillation in such a case reaches 
the maximum amplitude before a — 0, i.e, the red de- 
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tuning, corresponding to + Aq;2 < in Eq. (8). This 
implies negative coefRcients of quadratic and cubic terms 
in Eq. (4) originated from the different asymmetry from 
in z-axis. 



V. CONCLUSION 

In conclusion, we have experimentally investigated the 
complicated oscillations in our home-built SET, which 
are related to the high-order multipole potentials. Both 
the coupling and non-coupling cases, as well as the driv- 
ing along different axes, are studied. Our observation can 
be fully understood by the nonlinear effects and the mo- 
tional coupling from the solution of a complicated Duffing 
oscillator. 

In comparison to the relevant study on a single ion 
oscillating in the linear trap [isj . our home-made SET 
owns higher-order multipole potentials, which cause more 
fruitful nonlinear effects and even the motional coupling 
between different directions. Although there are also 
axial-radial couplings observed in the ion trap, e.g., with 
ion cloud in [2^ , such a motional coupling is much more 
evident in the SET, which, in addition to the couplings 
regarding zx and xz^ , is also reflected in the zx^ term in 
Eq. (4) with the coefficient a^. According to our calcula- 
tion, the motional coupling in our observation is mainly 
influenced by the coefficients as and ag, which implies 
the combined action from the quadrupole, hexapole and 
octopole potentials. This is the reason that the rectangle 



trajectories have never been observed previously in linear 
ion traps. Moreover, recent investigation of the phonon 
laser based on the nonlinear oscillation of the trapped 
ion demonstrated the analogy to the Fabry-Perot laser 
with 100% reflecting mirrors [l^. In contrast, the SET 
under our study seems an asymmetry Fabry-Perot cav- 
ity, which may yield two split beams of the phonon laser 
in perpendicular axes. Coherent transfer between the 
two split phonon beams would be useful in fundamental 
physics and practical application. Further work in this 
aspect is underway. 

With the trapped ions cooled down to the motional 
ground state, we may have an excellent platform to 
demonstrate nonlinear behavior in a fully quantum me- 
chanical regime and also carry out quantum logic gate op- 
erations. Therefore our work presents a way to exploring 
complicated nonlinearity using experimentally available 
techniques and it is also useful for suppressing detrimen- 
tal effects from the nonlinearity in quantum information 
processing using trapped ions. 
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Appendix I The nonlinear coefficients in Eq. (4) 

Our calculation is based on Eq. (4), in which the nonlinear coefficients originate from the high-order multipole 
potentials. By setting e and m to be the electric quantity and the mass of single calcium ion, we have the nonlinear 
coefficients as 

me^liUl^fMl^ ^ 140eM2i [mrtiV^*^^'' + 8er^ {Uj,) (t/f ) M,] ^^^^ 
6eroMi3 [mrliV*^)^-' + 6erg {U*) jU^,) Mj] ^^^^ 
3e [70erlMeM2i{U;)iU^,)+42erlMjM2oiUi)iUio)] 

OL2\ = 



a22 



^ 3e [2AerlMi2Mi3{U*^){U*^) + Imr^n^ M2o{V2*o)] 

3e [70erlMsM2i{Ui){Uii)+42er^oMYM22{U?){Ui2)] 
^ 3e [24ergMi3Mi4(t/r3)([/r4) + 7mr^nH422iV2*2)] 



(12) 



(13) 
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a4 



T* ^2 Ty f2 



e [-60mr^(F2*i)f^2M2i - Umr^iV^^)^^ M23 - 240er^{U^) ([/^i) M7M21 - 56er§(f77*) ([72*3) M7M23] 



,(14) 



0:5 



32er2 ([7*4)2 - 18er^ (JJl^y M'(^ + 6er^ ([/^g) ([/^s) M13M15 + 21eriMsM22iUI)iUi2) 



„2 /rr* \2 A,/2 



e [-60mr4(y2l)r!'A/2i + Um4{V2%)n''M23 - 240erg(£/f ) ([/2*i) AfrMzi + 56erg([/f ) ([72*3) Af7A/23] 



,(15) 



[6erg ([7*3) jU*^) Mi^Ahi + 14mrg(V,* )1]2A/i9 + 21erg ([/*) (L/2*2) M6M22] 

e [64erg([/i*2) (;7i*4) Afi2Afi4 + 56er§ (L/f ) ([/i*9) AfyMig + 2ler^„MsM2o{Ui){Uio)] 



(16) 



a7 



as 



[6erg ([/*3) ([/*) A/13A/6 + 8mrl{V{^)rtHh2 + 32erg ([7*2) jU^) A//12A/7] 
[8mr^(Vi*4)r!2A/,4 + S2erl ([7*4) (C/7*) A/i4Af7 + 6erlMsA'h3iU;)iU^3)] , 



with cjg^ = e 8er;5 ([7^)^ A/| + Amr^{Vj)n'^ Mr /m'^r^n'^, and the scahng factor tq (See definition in [22 



(17) 
(18) 



Appendix II Details of the steady-state solution Eq. (5) 



Eq. (5) is obtained by the standard steps of the multiple scale method. Star ting from Eq. (4), we assmne that the 
driving frequency is a perturbative expansion of harmonic oscillator frequency |27| . i.e., lUz = tooz + e^c, with a small 
dimensionless parameter e. Following the method in [28j . we rewrite Eq. (4) by setting z = eit, x = ep, y = eq, the 
damping term as 2c' ^ii and the driving term as t^k^ cos{uizt). Introducing a new parameter = eH (i — 0,1,2), we 
rewrite u,p, q as, 



u = e°uo {To, Ti,T2) + e^wi (Tq, Ti,T2) + £2^2 (Tq, Ti, T2) 
p = e% (To, Ti, T2) + eVi (To, Ti, T2) + e'p2 (Tq, Ti, r2) 
<Z = e"go (To, Ti, T2) + e^gi (To, Ti, T2) + ^q2 (To, Ti, r2) 

Then we compare the coefficients of e", and e^, which yields, 

T'o"o +1^02^0 = 0, 



(19) 



(20) 



Dlui + Wq^Wi = -2DoDiUa - a2ul - a-jq^UQ - asPoUo, 



(21) 



Dlu2 + ^02 ""2 



-[-kz cos (wo^To + CTr2) + qoa4Mo +Po"5ito +Po9oa6Wo + 2/iT'o"o 
+DIuo + 2D0D2U0 + goa2iw§ +poa22Wo + "3^0 + 2DqDiUi 
+0790^1 + asPoMi + 2a2MoUi + ar^iwo + asPiito], 



(22) 



where A = S/^T,, (i = 0, 1, 2). 

_We may solve uq = A (T2) exp(iwozTo) + A (r2) exp(-ia;o2To) and m = 02(^2 exp(2iwozTo) - QAA + 
exp(— 2ia;ozTo)]/3cjQ^ from Eqs. (20) and (21) by eliminating the secular term. Similarly, from equations of 
the oscillations in x-axis and y-axis, we may solve the variables p^, go, Pi and qi as 



Po = C (T2) exp(iiuoxTo) + C (T2) cxp(-iwoxTo), 



(23) 
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qo= B (T2) exp(ia;oj;To) + B (T2) cxp(-ia;oj;To), (24) 
Pi = a2AC^ exp(2ic^o.To) - 6CC + exp{-2iLUoccTo)]/{3Lul), (25) 

qi = a2y[B^ exp(2^woyTo) - 6BB + & exp(~2^L^oyTo)]/(3c^o\): (26) 

where a2x and a2y correspond to the quadric nonlinearity of the ion motion equation in x and y directions, respectively. 
a;oz:/27r and LOQy/2-K represent the harmonic frequencies in x direction and y direction. A = iaexp(i/3), B = i6exp(i<j) 
and C ~ ^cexp(z77), where a,b,c,/3,'^ and rj are real functions of T2, /3, <j and 77 represent the phases of different 
dimensions. A, B and C are conjugate terms of A, B and C. Substituting uq and ui into Eq. (22), we obtain an 
equation regarding the secular term, from which, in combination of Eqs. (23-26) with the expressions of A, B and C, 
we obtain 

1 . da 

- -kzSm'-f + afiujQ^ + ujo^—— = 0, (27) 



-k^ C0S7 -a4ab - -agac + — ^ + — + — + awo^ " l^^r = 0, (28 

2 8 4 4 12w;f^ 4a;^^^ iuj^y \ dT2 J 

with 7 = (TT2 — /3. We assume the steady-state motion corresponding to -^r ~ ~ 0. So we have 

f 3a3 3 3 , "4 , 2 "5 2 "7"2y ,2 "80:2^ 2 ^ , ^ \2 1 7,2 /on^ 

l~8~" "12^" '^T"^ "1^"*^ -aujozdj +{auJoz^l) ^ ^k^. (29) 

which is actually Eq. (5). 
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